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A relation giving a minimum for the irreversible work in quasi-equilibrium processes was 
derived by K. Sekimoto et al. (K. Sekimoto and S. Sasa, J. Phys. Soc. Jpn. 66 (1997), 
3326) in the framework of stochastic energetics. This relation can also be written as a type 
of "uncertainty principle" in such a way that the precise determination of the Helmholtz free 
energy through the observation of the work < > requires an indefinitely large experimental 
time At. In the present article, we extend this relation to the case of quasi-steady processes 
by using the concept of non-equilibrium Helmholtz free energy. We give a formulation of 
the second law for these processes that extends that presented by Sekimoto by a term of the 
first order in the inverse of the experimental time. We apply the results to a simple model. 

§1. Introduction 

In recent years, there has been a growing interest in studying small mesoscopic 
systems, immersed in different substrates, such as colloidal particles, nanoparticles 
in solutions, or biological systems, all of which are dominated by fluctuations. The 
principal interest is motivated due to recent experimental breakthroughs and techni- 
cal applications. Thermodynamic notions, such as applied work, dissipated heat and 
entropy, have been used successfully to analyze processes, in which single colloidal 
particles or biomolecules are manipulated externally.'^' In this area, it is possible 
to find several studies that have focused on the generation of non-equilibrium situ- 
ations, from a time-dependent potential, manipulated externally to model the effect 
of moving laser traps, micropippetes, or atomic force microscopic tips. In all these 
cases, it is straightforward to find a clear identification of external work, internal en- 
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ergy and dissipated heat, whose consistency has been proven in many experiments, 
going from the micro to the nano worldP''3'''21''Eli'El''El' On the other hand, systems 
that involve external hydrodynamic flow in soft matter systems have been studied 
theoretically using concepts like work, entropy production and with Langevin type 
equations of motionP' However, there is a deeper and more subtle issue hidden 
behind the proper treatment of external flows. Its presence or absence depends on 
the frame of reference. An adequate description of these systems should allow a 
frame- independent identification of the thermodynamic concepts 

The second law of thermodynamics describes the fundamental limitation on the 
possible transitions between equilibrium states. However, in contrast to equilibrium 
systems, the treatment and understanding of non-equilibrium systems requires going 
beyond the postulates of equilibrium thermodynamics. This unified framework al- 
lows the description of both the equilibrium and non-equilibrium phenomena, from 
an extension of the second law to a more inclusive state space of equilibrium and 
non-equilibrium steady states .'^''^''^ In this context, several authors have also 
established a connection between the phenomena related to non-equilibrium steady 
states and the thermodynamic laws for slow processes connecting different steady 

states.ram-ca 

The main objective of this paper is to shed some light on this last point, by 
addressing the analysis of a quasi-steady process, consisting of a particle in a moving 
substrate, while some control parameter changes slowly. A simple but paradigmatic 
case is a particle dragged through a viscous fluid by a harmonic laser trapP''^ with 
a quasi-statically changing strength. The dynamics of this particle can be modeled 
by a Langevin equation in the framework of Stochastic Energetics (see below). This 
approach will be done from the laboratory frame and the commoving frame, where 
the particle experiences a steady flow, and where a change of frame necessarily leaves 
the expressions for the thermodynamic magnitudes, as work and heat, invariant. 

According to thermodynamics, if we consider a system in contact with a heat 
bath and control parameters changes quasi-statically, the work W needed for the 
change is equal to the variation of Helmholtz free energy, AF: 

W = AF, (1-1) 

being AF composed by the sum of the reversible heat released to the heat bath, 
Qrev and the change of internal energy, AE: AF = Qrev + AE. If the change of the 
control parameters is not quasi-static, then the necessary work W is larger than the 
reversible one 

W - AF = Qirr > , (1-2) 

where Qirr is the irreversible heat that is equal to the difference between the 
released heat Q and the reversible heat: Qirr = Q — Qrev The released heat Q 
satisfies the first law : 



Q + AE = W. (1-3) 
In order to obtain the released heat Q for a given protocol of control parameters. 



Complementarity relation for quasi-steady states 



3 



it is necessary to have both a dynamical model of the system and a kinematical 
interpretation of the heat released by the system. An approach was introduced to 
obtain Q for systems whose dynamics are decribed by Langevin equations which 
are now known as stochastic energetics (SFi)^^^^ It constitutes an intermediary 
level of description that lies between Hamiltonian dynamics including all degrees of 
freedom of the concerned system, and thermodynamics where the system is controlled 
by external agents. In the framework of this approach, for a system that follows a 
quasi-static isothermal process, a complementarity relation giving a minimum for the 
product of the irreversible heat times the experimental time, Qirr^t > kBTSmin, was 
demontrated inP and an expression for the second law with a first order correction 
was obtained in]^ An extension of the stochastic energetics to the case of quasi- 
steady processes and an expression for the second law to zeroth order were presented 
jj^jigj jj^ ^j^g present letter we continue, in one sense, the work initiated irP^J by 
generalizing the approach used in^ We are able to show a complementarity relation 
that is valid for quasi-steady processes together with an expression for the second 
principle with a first order correction. This work is organized as follows: in the 
next section we apply the approach of SE to the case of the quasi-steady processes 
followed by a system satisfying a given Langevin equation (already used irP^ but 
slightly more general). At the same time, we sketch the principal steps of SE. We 
obtain our results and then exemplify them by studying a simple model. The next 
section is for the conclusions. Some computing is included in the appendix, in order 
not to deviate the text from the principal line of reasoning. 

§2. Stochastic energetics. The case of irreversible processes near steady 

states. The second principle 

We are going to extend the approach followed irP for irreversible processes 
near equilibrium states, to the case of irreversible processes near steady states. We 
study the response of a system in contact with a single bath, at temperature T, 
in the limit of very strong friction. Let x = {xi,...Xn} representing the state of 
the fluctuating system and b = {bi,...br} the parameters that control the system 
through the potential [/(x(t); a(t); b(t)) = t/(x(t) - a(t);b(t)). The quantity a{t) 
is another parameter that models a conservative force, depending on the distance 
to a major object, or bound to an arbitrary origin, becoming relevant if we make a 
Galilean transformation, as we will see later. Besides the conservative forces arising 
from the potential C/(x — a), we can perturb the particle by a direct force f. This 
force may include all the conservative and non-conservative contributions that do 
not depend on a. We describe the stochastic particle motion through the Langevin 
equation as follows: 



where is a friction constant given by a symmetric and positive definite matrix, 
and ^{t) is a Gaussian and white-correlated stochastic force, satisfying 
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where P = t^. 

In close analogy with the case of quasi-equilibrium process studied infD we 
rewrite Eq. ()2-ip by making the scalar produclCJ by dx along the realized trajec- 
tory and using that dU = ^-dx + + ^da, obtaining a balance equation for 
energy, which is: 



^ f/x \ 8U 8U 

r- — -${t)-f \ •(ix + dC/(x;a;b) = — (x;a;b)-db + — (x;a;b)-da (2-3) 

or re-ordering terms: 



fr— -^(t)Vdx+dC/(x;a;b) = ^ (x; a; b) •db+?^ (x; a; b)-da+f -dx . (2-4) 
\ at J OD oa 

We can rewrite this last expression as 

dQ + dU = dW , (2-5) 

where 

^iQ= (r.^-^(t)) -dx (2-6) 
is the heat discharged onto the bath and 

dW = ^ (x; a; b) • db + ^ (x; a; b) • da + f • dx (2-7) 
oh oa 

is the total work done by the external agent to the system. 

It is possible to visualize the problem from another reference frame, S'{X,t), 
which is obtained by making a Galilean transformation from x to X, taking a{t) = 
v.t, by mean of 

X = X-Fa(t) = X-h v.t (2-8) 

which implies 

± = X + v, (2-9) 
where '= ^. Transforming the Langevin equation ()2-ip in S'{X,t), we obtain 

^ ■ ? = - + ^(*) + f (2-10) 

where the potential in the new reference frame /i(X;b), depending now on X 
and on the parameters h h : i — > R, is obtained from the potential f7(x; a; b),as: 



*-* The multiplication of fluctuating quantities, i.e. ^{t)) ■ dx, should be understood in the sense 
of Stratonovich calculus.'^ 
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(x; a; b) A (x — a; b) A [/(x; a; b), i.e.: 

/i(X; b) = U{g{^; a; b)) = [/(x - a; b) (2-11) 
where g : R"^ i — > is defined as g{^; a; b) = (x — a; b) . 

Thus, we study the response of the system over the realized trajectory with 
respect to a steady state. This state will be reached by the system if it is not 
perturbed and it is uniquely determined by the parameter values (i.e. a) 

It is easy to verify that: 



dh _ dU 
9X ~ 9l ' 

dh _ dU 
5b ~ 9b ' 

dh _ dU_ 
5X ~ ~ 9a ■ 

Thus, for the potential variation we verify the equalitj0 



(2-12) 
(2-13) 
(2-14) 



dU = dh . (2-16) 
The energy balance equation reads 

r • ^ - ■dx + dh = |^(x(i); b) • db + (/ - rv) ■ dx (2-17) 

that we write as 

dQ' + dU = dW' , (2-18) 
where we have defined the heat and the work in the S' system as: 

dQ' = ir-^-m)-dX (2-19) 

dh 

dW = ^(X(t); b(t))db + (/ - rv)dX . (2-20) 

Equations (j2-5p and (|2-18p constitute the expression for the energy balance law 
in two different reference frames, related by the Galilean transformation (|2-8|) . We 
see that the law kept its same invariant form under a Galilean transformation, in 
agreement with the first law: dU = dQ — dW. It is easy to verify that the laws of 
transformation for the heat and for the work are given by: 



*■* This can be deduced from the invariance under Gahlean transformations of the scalar product 
dU . dU dU . dh .„ dh 
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dh 

dQ = dQ' + {-— + f -rv)dai + rvdjc , (2-21) 
oX. 

dh 

dW = dW' + + f - rv)da. + Evdx . (2-22) 

Taking the average of the quantities in Eqs. ()2-2ip and ()2-22p in the steday state, 
that is considering the probabihty distribution function of the steady statCj for a 
given dt, we get 

< dQ >=< dQ' > +rvdj^ (2-23) 

< dW >=< dW' > +rt;(ix , (2-24) 
where the identity < (— ^ + / — -^v) >steady = was used. 
The term Fvdx represents the housekeeping work. 

The work {W) in the frame S'(X,t). 

If the control parameters b change from b(0) = bj to h{At) = bj, then the total 
work performed on the system along a particular process X(t) (0 < i < At) is given 
by 

W = j\t^{y.{t)-Mt)) • ^ + |</>dx , (2-25) 

where we have defined (p = f — Fv. 

The ensemble average of the work, (W), over a possible realization of {^(i)}o<t<z\t 
can be computed as 



(W) - { I (/.dX) = l\t I dXP(X,t)|^(X;b(t)) 



where P is the probability distribution function of X that satisfies the Fokker-Planck 
equation, which is given by 

^(X, t) = -CFPih{t))PiX, t) (2-27) 

where 



CFP{h{t)) ^ A . • * (^(X; b(t)) + rv(t) - f + ksT^^ . (2. 



28) 



*■* The probability distribution is given in Eq. (|A-7|I . 
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<W'>-< [(l)dX>= AF*+ / < X > d<P+^ C ds^^yl(b)^^ , At ^ oo , 
J J At Jq ds ds 

(2-29) 

where F* is the non-equihbrium free energy, defined in Eq. (|A-13P and ^(b) is 
a positive definite n x n matrix defined in the appendix, Eq. ()A-18p . 

By means of a Legendre transformation, we define the new free energy G* (T, (X) , b) 

as 

G*(T, (X),b) ^ F*(T,0,b) - ^^*^gf'^^ ■ <P (2-30) 

with 

dF*(T,6,h) 

\ V^' ^ = - <X > 2-31 

Then we have from Eq. (|2-29|) 

<W' >+{[ <pd<X>-< UdK>) = AG*+^ C ds'^^A{h)^^ , At ^ oo , 
J J At Jq ds ds 

(2-32) 

and using that < (j) >= (j) 

<W'> + < /" 0(d < X > -dX) >= AG*+^ C ds^^^A{h)'^^ , At ^ oo . 
J At Jq ds ds 

(2-33) 



The term < J < X > — dX) > is equal to zero to &^orde7**\ 

Then the total irreversible work, iyV') — AG* , for a very slow process is given 

by 

The integral on the r.h.s of Eq. ()2-34p has the form of a classical action for a 
particle of "mass" /1(b) and has a minimum Sxnin{,^i^^ for a certain "classical 
path. Hence, as in the case of quasi-equilibrium process, an inequality that resembles 
a sort of " uncertainty" relation remains true for the present case of steady process, 
valid for Z\t — )■ oo: 

{{W) - AG*) At>Sm^n{h^,hf) . (2-35) 

According to ()2-35p , the estimation of the non-equilibrium Helmholtz free energy, 
by the measurement of the net mean mechanical work, contains an indetermination 



*-* See the appendix where, in order not to deviate the text from the principal line of reasoning, 
the computations are provided. 

**' < d(< X > —X) in each instant if the parameter b changes very slowly(0*'' order) 
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Qirr = {W) — AG* (the total irreversible work), whose product by At cannot 
be smaller than a positive lower bound. The precise determination of the non- 
equilibrium Helmholtz free energy through the observation of the work (W) requires 
an indefinitely large experimental time At. 

We can express our results in differential form, i.e., for an elementary process. 
From (|2-34p . we have up to the first order 

{dW) = dG* + ^A{h)^dt (2-36) 

Lit/ (JiL 

where we used (|A-1|) to return to variable t. Equation ()2-36p represents the 2nd law 
for quasi-steady processes with a 1st order correction. We were able to obtain a 1st 
order correction to equation (5-3), which corresponds to the case of a simple 

model for the steady-state thermodynamics presented there. 



Application to a simple model. 

As an application of the preceding approach for steady processes, we consider a 
single particle trapped in an one-dimensional harmonic potential U{x, b{t)) = ^6(t)x^ 
being the strength b{t) the control parameter that changes slowly during a time 
lapse. The variable x represents the displacement of the particle from the origin. It 
is immersed in viscous external flow with velocity v at temperature T; F denotes 
the friction constant. We want to calculate the total work done to the system, in the 
course of the process, concerning the change of b, from different reference inertial 
systems S (laboratory frame) and «S" (comoving flow frame). The relevant parameter 
for us is the flux velocity v, or the relative velocity, between the two references frames. 
Thus, we can redefine the potential in S' as h{X{t),b{t)) = ^b{t)X^. 

In the absence of an external force (f = 0), the differential expression for the 
excess of work dW obtained from Eqs. ()2-2ip . ()2-22p . and the corresponding Langevin 
equation takes the following form: 

{dW) = {dW) - ^ (^^"j'dt (2-37) 
that, according to Eq. ()2-36p is equal to: 

(dW) = dG* + ^A{b)^dt . (2-38) 

This expression is the 2nd law for the quasi-steady isothermal processes, with a 
1st order correction for this paradigmatic model. If we discharge the first order term 
in the r.h.s., the last expression reduces to Eq. (5-3) of Sekimoto.'^ Furthermore, 
we see that the housekeeping work appears to be naturally subtracted from the mean 
work. 

We were able to compute the quantity ^(b) from its simplified formula Eq 
(|AT9p . now that (one-dimension) is a 1 x 1 matrix, i.e. a scalar quantity. In this 
model, the steady distribution is given by 
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-I3i^b(t)x''+rvx) 



and the kernel g(X., X'; b) is 



ff(X, X'; b) = ^sgn{X - X') T ^ . (2-40) 



2 -^-'-^^ ^^'Jx, Pst{X;b,v 
After integration of ()A-19p . we found: 



r r^v 



2 



thus, from ()2-38p . we have for the irreversible heat: 



dQirr = ( T^T + ^ ) b^dt (2-42) 



^4/363 64 

We see that two quite different terms contribute to A{b), and therefore to the 
irreversible heat released during the whole process of variation of the parameter b. 
By means of qualitative analysis, we are going to show that the difference between 
these two terms lies in their physical origin. The first term depends directly on 
temperature, and contains the information about the effect of the change of b on the 
mean fluctuations of the particle position, when it is coupled to the thermal bath. 
The second term is the only one with a contribution due to the steady regime, that 
is, depends on the velocity v, and is due to the rearrangement of the equilibrium 
position Xeq. 
In thermal equilibrium, 

hsT =< h{x) >= < X2 > , (2-43) 

thus, we can define a typical mean displacement (or equivalently, the typical 
mean amplitude of oscillation) as: 



= (2-44) 



with an equilibrium position that in the steady regime is 

X,g = -^ . (2-45) 

We associate Xa and X^q with typical velocities of adjustment, vi V2, due to the 
change of the parameter b. They are given by 

vi=Xd = -lVk^^ (2-46) 



2' 65 
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and 

Ev ■ 

V2 = = -^h . (2-47) 

Then we can estimate the irreversible heat as the energy dissipated by the friction 
force, Evi, (i = 1,2), in the viscous medium, in the infinitesimal lapse dt. It is given 
by 

dQ^rr = Evldt + Evldt = f ^ + j h^dt (2-48) 

which is exactly Eq. (|2-42p . 

It is important to note that for the case f = (quasi-equilibrium) we re-obtain 
the result found for A in Eq.(18) of RefP When A is replaced in expression Eq. 
(j2-36p and is integrated during the time lapse, the integral can be minimalized in 
order to obtain an optimal protocol for which the irreversible heat is minimumP 

§3. Conclusion 

Following an analogue approach for the case of the quasi-equilibrium processes 
we were able so show that an inequality, connecting the irreversible work (W') — AG* 
and the experimental time At, exists for the case of quasi-steady processes. It is given 
by ()2-35p and it states that the estimation of the non-equilibrium Helmholtz free 
energy, by the measurement of the net mechanical work, contains an indetermination 
Qirr = (^') ~ AG* (the total irreversible work), whose product by At cannot 
be smaller than a positive lower bound. The precise determination of the non- 
equilibrium Helmholtz free energy through the observation of the work {W') requires 
an indefinitely large experimental time At. We deduced the second law for quasi- 
steady processes with a 1st order correction for a particle subject to a potential 
and immersed in a flow with velocity and we showed that the HKW appears to be 
naturally subtracting from the mean work. 
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Appendix A 

Proof of Eq. 12- 29\) : the work in the moving frame. 

In order to compute the integral from (|2-26p . the scaled time s 

*-* The present potential was studied inj^ Remark 4, where the elastic constant is named a 
instead of b. 
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At 



(A-l) 



as 



is defined. The probability distribution depending on this argument is defined 
P(X,s;Z\t) = P{X,sAt), 

and the parameters as b(s) = h(sAt). Equations ()2-26p and ()2-27p become 



(W) 



(/.dX) 



'ds^. /dx|^(X;b(5))P(X,s;Z\t), 







ds 



dh 



l-^{X,s;At) = -CFp{Hs))P(^,s;At) . 



(A-2) 



(A-3) 



Eq. ()A-3p can be solved perturbatively by assuming that At is large enough to 
make an expansion of P in powers of as 



P(X, s;At) = (X, s) + i^p(i) (X, s) + • • • 



Substituting in ()A-3p . we have for the zero and first order 

= -£Fp{h{s))P^°\x,s), O^'^order 



ds 



(X, s) = -CFp{h{s))P^^'^ (X, s) reorder. 



(A-4) 



(A-5) 



(A-6) 



From the lowest order, Eq. (|A-5p . and the normalization condition J dXP^'') (X, s) 
1 we deduce that P^'^^ is the steady distribution Pst for a given parameter b(s): 



P(°)(X,.) = P.,(X;b(.),v) ^ ^^^^_,(,(x,b)4-^.x) 



(A-7) 



,Kp and (j) = f — Fv. 
If f = cte, Eq. ()A-6p becomes 



where /3 = 



dP. 



St 



ds 



:X;b(s),v) = -/:^p(b(s))P«(X,s) 



Now, the kernel ^(X, X'; b(s)) is defined as the solution of 

-CFp{h) rP,t(X;b(s))<7(X,X';b(s))l =5(X,X') 



(A- 



(A-f 



If we multiply Eq. ()A-8P by P^f (X; b(s))(7(X, X'; b(s), v) and then integrate in 
X, we obtain p(^)(X, s) as 



P«(X,s) =P,t(as)(X;b(s),v) 



dX'5(X,X';b(s)) 



^(X';b(s),v)+x 



(A-IO) 
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where the integration constant x is obtained from the normahzation condition, 
/ ciXP(i)(X,s) =0, as 



X 



dx<|P,i(X;b(s)) / dX'5(X,X';b(s))^(X';b(s)) 



(A-ll) 



Having obtained -P(X, s) up to the first order, we substitute 

P(X,s) = P,t(X;b(s),X) + ^p(i)(X,s) + ... in Eq. and we have 

(A-12) 

As we are deahng with an out of equihbrium process, it is useful to make use 
of the non-equihbrium Helmholtz free energy F*{T, b, v), defined by Sekimoto in^^ 
that is given by 



F*(r,0,b) = -ksTln 



exp - 



/i(X;b) 



•X 



-dX 



(A- 13) 



where (/> = f — ^^^ and /i(X; b) — (/> • X is an "effective" potential. 


The following "Ehrenfest type" identity, concerning the steady ensemble average 
(Ie>p,,' is satisfied: 

dF* 



dh\ 



and furthermore 



dF* 



dh - 
dX— (X;b)P,t(X;b) , 



<X> , 



(A-M) 
(A-15) 



so we have, from Eq. (jA-12p . up to the first order 



{W')-{ <pdX) 



(A-16) 

Using the relation ("'chains rule"') ^(X';b,v) = * (^^(X'; b(s), v)) • 

and substituting ()A-10p (using ()A-lip ) in the first order term of ()A-16p . we have for 
{W) 



{W')-{J </'dX) = ^ (^dh 



^ f „ dF* 1 /"\ (ib(s) ,db(s) , 2n 



dh 



At 



ds 



ds 



(A-17) 



*■* We can obtain F*(T, 0, b) in an operational way by mean of the equation (|2-29|) . 
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where 



.1(b) ^ j dxj dX!' 

j dx(5(X-X)-P,i(X;b,v))5(X,X';b)(^^(X';b,v)) (A-IS) 



is a positive definite n x n matriJ^. which can be simpHfied to 



Aih) = -^Jdxl dX' ^(X;b,v).5(X,X';b) * (^^(X';b,v)) .(A-IO) 



being 5(X,X';b) the Green's function satisfying Eq. (|A-9p . 

From the definition of the non-equihbrium Helmholtz free energy, Eq. (jA-lSp . 
for T = constant, we have 

dF* dF* dF* 

Substituting ^AJW^ in (lA47j) it follows Ea. (l2^ . 
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